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Abstract 

We review and update the constraints on the parameters of the quark flavour mixing 
matrix Vckm in the standard model and estimate the resulting CP asymmetries in B 
decays, taking into account recent experimental and theoretical developments. In per- 
forming our fits, we use inputs from the measurements of the following quantities: (i) |e|, 
the CP-violating parameter in K decays, (ii) AM^, the mass difference due to the B®-B® 
mixing, (iii) the matrix elements \V c b\ and \V u b\, (iv) 5-hadron lifetimes, and (v) the top 
quark mass. The experimental input in points (ii) - (v) has improved compared to our 
previous fits. With the updated CKM matrix we present the currently-allowed range of 
the ratios |Vtd/Vt s | and |Vtd/^u&|, as well as the standard model predictions for the B®-B® 
mixing parameter x s (or, equivalently, AM S ) and the quantities sin 2a, sin 2/? and sin 2 7, 
which characterize the CP-asymmetries in I?-decays. Various theoretical issues related to 
the so-called "penguin-pollution," are of importance for the determination of the phases 
a and 7 from the CP-asymmetries in B decays, are also discussed. 



*Presented at the 6th. International Symposium on Heavy Flavour Physics, Pisa, June 6 - 10, 1995. 



1 Introduction 



The aim of this article is to revise and update the profile of the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix reported earlier by us in particular the CKM unitarity tri- 
angle and the CP asymmetries in B decays, which are the principal objects of interest 
in experiments at present and forthcoming B facilities. In performing this update, we 
include the improvements reported in a number of measurements of the lifetime, mixing 
ratio, and the CKM matrix elements \Va\ and \ V u t,/V c b\ from B decays, as well as the top 
quark mass. On the theoretical side, we mention the improved estimates of the power 
corrections in the analysis of the exclusive semileptonic decay B — > D*lvi in the context 
of the heavy quark effective theory (HQET) || [|, and the calculation of the missing part 
of the next-to-leading order calculations in the analysis of the CP- violating quantity |e| 
0]. We note here the changes that we have made in the input to our present analysis 
compared to that reported by us in Ref. 

• The top quark (pole) mass m t = 174 ±16 GeV, measured earlier by the CDF collab- 
oration |J, is now replaced by improved measurements by the same collaboration 
|| and by DO @, yielding the present world average m t = 180 ±11 GeV 0. This 
leads to the running top quark mass in the MS scheme, m7(m t ) = 170 ± 11 GeV 



A new and improved measurement of the quantity jF(l)|V^h| in the decays B — > 
D*tv^ using methods based on the heavy quark effective theory (HQET), has been 
reported by the ALEPH collaboration |fL0[j . This is lower than their previous number 



]nj] , as well as the corresponding numbers from the CLEO [|T2j and ARGUS [|Lj 
analyses. Likewise, new measurements are reported by the DELPHI collaboration 
[[bill . In the meantime, estimates for the quantity T{V) = £(1)77,4 have undergone 
some revision in both the QCD perturbative part tja and power corrections to the 
Isgur-Wise function at the symmetry point £(1). We use the value ^(l) = 0.91 ± 
0.04, obtained recently by Neubert [Q, and which is in good agreement with the 
estimates of Shifman et al. || Taking into account the updated experimental 
and theoretical input, we obtain \V c b\ = 0.0388 ± 0.0036. The central value for this 
matrix element has come down compared to the value \V c b\ = 0.041 ± 0.006 used by 
us previously, and the error on this quantity is now smaller, about ±9%. 

Until recently, the knowledge of the CKM matrix element ratio [V^/V^l was based 
on the analysis of the end-point lepton energy spectrum in semileptonic B decays 



nj, which is quite mo del- dependent. We had used a value [l^b/K&l = 0.08 ± 0.03 
to take this model dependence into account. In the meantime, the measurement of 
the exclusive semileptonic decays B — > (ir, p)tvt has been reported by the CLEO 
collaboration [17]. The matrix element ratio so determined is also model-dependent 
due to the decay form factors. However, this set of data permits a discrimination 
among a number of models, all of which were previously allowed from the inclusive 
decay analysis. The convolution of the two methods reduces the theoretical uncer- 
tainty somewhat. We use a value |T4&/K:&| = 0.08 ± 0.02, which is a fair reflection 
of the underlying present theoretical dispersion on this ratio. 
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In the analysis of the CP- violating quantity |e|, the perturbative renormalizations of 
the various pieces in the | AS\ = 2 Hamiltonian from the intermediate charm and top 
quark are required |18| . While the next-to-leading order results for the quantities 
f] cc and fju have been known for some time and were used in our previous analysis, 
the next-to-leading order calculation of the quantity f] ct has been completed only 
recently Q. We use the improved calculation of r) ct in performing the CKM fits 
presented here. 

The measurements of the B®-B® mass difference AM^ have become quite precise, 
with the present world average being AMj = 0.465 ±0.024 (ps)^ 1 [[19]. The present 



lower limit on the B° s -B® mass difference has improved slightly, AM S > 6.1 (ps) 1 
at 95 % C.L., assuming for the probability of the fragmentation of a b quark into a 
B s meson a value f s = 12% H, yielding AM s /AM d > 12.3 at 95% C.L. gj|. 



All of these improvements warrant an updated fit of the CKM parameters. 

As in our previous analysis, we consider two types of fits. In Fit 1, we assume particular 
fixed values for the theoretical hadronic quantities. The allowed ranges for the CKM 
parameters are derived from the (Gaussian) errors on experimental measurements only. 
In Fit 2, we assign a central value plus an error (treated as Gaussian) to the theoretical 
quantities. In the resulting fits, we combine the experimental and theoretical errors in 
quadrature. For both fits we calculate the allowed region in CKM parameter space at 95% 
C.L. We also estimate the SM prediction for the B° s -B® mixing parameter, x s , and show 
how the ALEPH limit of AM s /AM d > 12.3 (95% C.L.) constrains the CKM parameter 
space. We give the present 95% C.L. upper and lower bounds on the matrix element ratio 
|Vtd/Vf S |, as well as the allowed (correlated) values of the CKM matrix elements \V t d\ and 
\V ub \. 

We also present the corresponding allowed ranges for the CP-violating phases that 
will be measured in B decays, characterized by sin 2(3, sin 2a and sin 2 7. These can be 
measured directly through rate asymmetries in the decays Bd — > J/ipKs, Bd — > 7r + 7r~, 
and B s — > DfK T (or B — > D K ), respectively. We also give the allowed domains for 
two of the angles, (sin 2a, sin 2(3). Finally, we briefly discuss the role of penguins (strong 
and electroweak) in extracting the CP phases a and 7 from the measurements of various 
CP-asymmetries. 

This paper is organized as follows. In Section 2, we present our update of the CKM 
matrix, concentrating especially on the matrix element \V c b\ which, thanks to the progress 
in HQET and experiments, is now well under control. The constraints that follow from 
l^u&Ac&l, |e| and AMd on the CKM parameters are also discussed here. Section 3 contains 
the results of our fits. These results are summarized in terms of the allowed domains of 
the unitarity triangle, which are displayed in several figures and tables. In Section 4, 
we discuss the impact of the recent lower limit on the ratio AM s /AMd reported by the 
ALEPH collaboration on the CKM parameters and estimate the expected range of the 
mixing ratio x s in the SM based on our fits. Here we also present the allowed 95% C.L. 
range for | V*rf / V* s | . In Section 5 we discuss the predictions for the CP asymmetries in the 
neutral B meson sector and calculate the correlations for the CP violating asymmetries 
proportional to sin 2a, sin 2(3 and sin 2 7. We also review some of the possible theoretical 
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uncertainties in extracting these CP-phases due to the presence of strong and/or elec- 
troweak penguins. We present here the allowed values of the CKM matrix elements \Vtd\ 
and \V u b\. Section 6 contains a summary and an outlook for improving the profile of the 
CKM unitarity triangle. 



2 An Update of the CKM Matrix 

In updating the CKM matrix elements, we make use of the Wolfenstein parametrization 
|2l| , which follows from the observation that the elements of this matrix exhibit a hierarchy 
in terms of A, the Cabibbo angle. In this parametrization the CKM matrix can be written 
approximately as 



(1) 



'CKM 



t 1 - |A 2 

-A(l +iA 2 X 4 v ) 1 
\AA 3 (1 -p-irj) - 




In this section we shall discuss those quantities which constrain these CKM parameters, 
pointing out the significant changes in the determination of A, A, p and r\. 

We recall that \V US \ has been extracted with good accuracy from K — > -aev and hyperon 



decays |22| to be 

(2) \V US \ = A = 0.2205 ±0.0018 . 

This agrees quite well with the determination of V u d ~ 1 — \\ 2 from /?-decay, 

(3) \V ud \ = 0.9744 ±0.0010 . 

The parameter A is related to the CKM matrix element V c b, which can be obtained 
from semileptonic decays of B mesons. We shall restrict ourselves to the methods based 
on HQET to calculate the exclusive and inclusive semileptonic decay rates. In the heavy 
quark limit it has been observed that all hadronic form factors in the semileptonic decays 
B — > (D, D*)£i>£ can be expressed in terms of a single function, the Isgur-Wise function 
| |23| |. It has been shown that the HQET-based method works best for B — > D*lv decays, 
since these are unaffected by l/mg corrections [0, ^5], |26|. This method has been used 



by the ALEPH, ARGUS, CLEO and DELPHI collaborations to determine £(l)\V cb \ and 
the slope of the Isgur-Wise function. 

Using HQET, the differential decay rate in B — > D*lv^ is 



dT{B -> D*tv) G 



2 



(4) ^ - ^ = ^{m B - mn.fmWA^^^ + if 



, 4a; 1 — 2ujr ±r 2 llT _ 
X|1 + ^TT (l-r)» 



where r = m^* /ms, to = v ■ v' (v and v' are the four- velocities of the B and D* meson, 
respectively), and t]a_ is the short- distance correction to the axial vector form factor. In 
the leading logarithmic approximation, this was calculated by Shifman and Voloshin some 
time ago - the so-called hybrid anomalous dimension p7| . In the absence of any power 
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corrections, £(u = 1) = 1. The size of the 0(1 /ml) and 0(1 /ml) corrections to the Isgur- 
Wise function, £(o>), and partial next-to-leading order corrections to t}a have received a 
great deal of theoretical attention, and the state of the art has been summarized recently 
by Neubert and Shifman ||. Following them, we take: 

= 1 + 5(l/m 2 ) = 0.945 ±0.025 , 

(5) T] A = 0.965 ± 0.020 . 

This gives the range [[|: 

(6) T(l) = 0.91 ±0.04 . 

The present experimental input from the exclusive semileptonic channels is based on the 
data by CLEO |[2j, ALEPH JTTJ, ARGUS @, and DELPHI 0: 



iVrtl-^l) = 0.0351 ±0.0019 ±0.0020 [CLEO], 

= 0.0314 ±0.0023 ±0.0025 [ALEPH], 

= 0.0388 ±0.0043 ±0.0025 [ARGUS], 

(7) = 0.0374 ±0.0021 ±0.0034 [DELPHI], 

where the first error is statistical and the second systematic. The ARGUS number has 
been updated by taking into account the updated lifetimes for the B° and B ± mesons 
The statistically weighted average of these numbers is: 



(8) \Vcb\ -F(l) = 0.0353 ±0.0018 , 
which, using T(l) from Eq. (0), gives the following value: 

(9) \V cb \ = 0.0388 ± 0.0019 (expt) ± 0.0017 (th). 

Combining the errors linearly gives |V^,| = 0.0388 ± 0.0036. This is in good agreement 
with the value \V cb \ = 0.037lg;^ obtained from the exclusive decay B — > D*£ue, using a 
dispersion relation approach [E?J. Likewise, the value of \V c b\ obtained from the inclusive 
semileptonic B decays using HQET is quite compatible with the above determination 



(10) \V cb \ = 0.040 ± 0.0010 (expt) ± 0.005 (th) . 
In the fits below we shall use \V cb \ = 0.0388 ± 0.0036, yielding 

(11) A = 0.80 ± 0.075 . 

The other two CKM parameters p and rj are constrained by the measurements of 
iKifi/Kbl) N (the CP-violating parameter in the kaon system), x rf (B^-B^ mixing) and 
(in principle) e'/e (AS = 1 CP- violation in the kaon system). We shall not discuss 
the constraints from e'/e, due to the various experimental and theoretical uncertainties 
surrounding it at present, but take up the rest in turn and present fits in which the allowed 
region of p and rj is shown. 

Up to now, \V ub /V cb \ was obtained by looking at the endpoint of the inclusive lepton 
spectrum in semileptonic B decays. Unfortunately, there still exists quite a bit of model 
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dependence in the interpretation of the inclusive data by themselves. As mentioned earlier, 
a recent new input to this quantity is provided by the measurements of the exclusive 
semileptonic decays B — > (ir, p)£u£. The ratios of the exclusive semileptonic branching 
ratios provide some discrimination among the various models |I7| |. In particular, models 



such as that of Isgur et al. pOU , which give values in excess of 3 for the ratio of the 
decay widths r(_B° — > p~£ + u)/T(B° — > 7r~£ + z/), are disfavoured by the CLEO data. The 
disfavoured models are also those which introduce a larger theoretical dispersion in the 
interpretation of the inclusive B — > X u £u£ and exclusive decay data in terms of the ratio 
l^u&Ac&l- Excluding them from further consideration, measurements in both the inclusive 
and exclusive modes are compatible with 

(12) — = 0.08 ±0.02 . 

V c b 

This gives 

(13) ijp 2 + r] 2 = 0.36 ± 0.08 . 

With the measurements of the form factors in semileptonic decays B — > (tt, p, lS)Iv^ one 
should be able to further constrain the models, thereby reducing the present theoretical 
uncertainty on this quantity. 

The experimental value of |e| is |22| 
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(14) |e| = (2.26 ±0.02) x 10 

Theoretically, |e| is essentially proportional to the imaginary part of the box diagram for 
K°-K° mixing and is given by ]3TJ 



(15) + %mh(v,)A 2 \ l (i-P)), 

where yi = m^/M^, and the functions f 2 and f 3 can be found in Ref. Here, the rji 



are QCD correction factors, of which f] cc [[32fl and f] tt |33[ were calculated some time ago 
to next-to- leading order, and f\ ct was known only to leading order |TS[ Q . Recently, this 
last renormalization constant was also calculated to next-to-leading order |4| . We use the 
following values for the renormalization-scale-invariant coefficients: f] cc ~ 1.32, fj tt ~ 0.57, 
fj ct ~ 0.47, calculated for rh c = 1.3 GeV and the NLO QCD parameter Ajjg = 310 MeV 
in Ref. §. 

The final parameter in the expression for |e| is the renormalization-scale indepen- 
dent parameter Bx, which represents our ignorance of the hadronic matrix element 
(iT°|((i7 M (l — 75)5) \K°). The evaluation of this matrix element has been the subject 



of much work. The earlier results are summarized in Ref. |BS| 



In our first set of fits, we consider specific values in the range 0.4 to 1.0 for Bk- As we 
shall see, for Bk = 0.4 a very poor fit to the data is obtained, so that such small values 
are quite disfavoured. In Fit 2, we assign a central value plus an error to Bk- As in our 
previous analysis [|l], we consider two ranges for B K - 

(16) £?^ = 0.8±0.2, 
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which reflects the estimates of this quantity in lattice QCD |31| [37], or 
(17) B K = 0.6 ± 0.2 , 

which overlaps with the values suggested by chiral perturbation theory p8[ . As we will 
see, there is not an enormous difference in the results for the two ranges. 

We now turn to B° d -B Q d mixing. The present world average of Xd = AMd/Td, which is 



a measure of this mixing, is [19 



x d = 0.73 ±0.04 



which is based on time-integrated measurements which directly measure Xd, and on time- 
dependent measurements which measure the mass difference AMd directly. This is then 
converted to Xd using the B d lifetime, which is known very precisely (r(Bd) = 1.57 ± 
0.05 ps). From a theoretical point of view it is better to use the mass difference AMd, as 
it liberates one from the errors on the lifetime measurement. In fact, the present precision 
on AMd, pioneered by time-dependent techniques at LEP, is quite competitive with the 
precision on Xd- The LEP average for AMd has been combined with that derived from 



time-integrated measurements yielding the present world average |19 
(19) AM d = 0.465 ± 0.024 (ps)' 1 . 



We shall use this number instead of Xd, which has been the usual practice to date |L8|, [35 



38, 39 



The mass difference AMd is calculated from the B d -B d box diagram. Unlike the 
kaon system, where the contributions of both the c- and the t-quarks in the loop were 
important, this diagram is dominated by t-quark exchange: 

(20) AM d = ^M 2 W M B (f 2 Bd B Bd ) f, B y t f 2 (y t )\V td V tb \ 2 , 

where, using Eq. [j], |V^V^| 2 = A 2 X 6 (1 — p) 2 + rf . Here, f] B is the QCD correction. 



In Ref. |33[, this correction is analyzed including the effects of a heavy t-quark. It is 
found that f) B depends sensitively on the definition of the t-quark mass, and that, strictly 
speaking, only the product fi B (y t )f2(yt) is free of this dependence. In the fits presented 



here we use the value f] B = 0.55, calculated in the MS scheme, following Ref. |33| . 



Consistency requires that the top quark mass be rescaled from its pole (mass) value of 
m t = 180 ± 11 GeV to the value rfTi(m t {pole)) in the MS scheme, which is typically about 
10 GeV smaller §. 

For the B system, the hadronic uncertainty is given by f B B Bd , analogous to Bk in 
the kaon system, except that in this case, also f Bd is not measured. In our fits, we will 
take ranges for f Bd B Bd and B Bd which are compatible with results from both lattice-QCD 
and QCD sum rules ft|, g0|, |4|: 



f Bd = 180 ± 50 MeV , 
(21) B Bd = 1.0 ±0.2. 

In Table [I], we summarize all input quantities to our fits, of which seven quantities (\V c b\, 
\V u b/V c b\, AMd, i~(Bd), mi, fj cc , fj ct ) have changed compared to their values used in our 
previous fit [jlj. 
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Pq vq ryipf pT* 
J. OjX ( X.lLlvZ tCI 


Vfl 1 n p 


A 


0.2205 




0.0388 ± 0.0036 




0.08 ± 0.02 


e 


(2.26 ± 0.02) x 10 6 




(0.465 ± 0.024) {ps)~ l 




(1.57 ± 0.05) (ps) 


m t [m t \poie)) 


(170 ±11) GeV 


Vb 


0.55 


Vcc 


l.oz 


'let 


0.47 


Vtt 


0.57 


Bk 


0.8 ±0.2 


Bb 


1.0 ±0.2 


fs d 


180 ± 50 MeV 



Table 1: Parameters used in the CKM fits. Values of the hadronic quantities fs d , Bb a 
and Bk shown are motivated by the lattice QCD results. In Fit 1, specific values of these 
hadronic quantities are chosen, while in Fit 2, they are allowed to vary over the given 
ranges. (In Fit 2, for comparison we also consider the range Bk = 0.6 ± 0.2, which is 
motivated by chiral perturbation theory and QCD sum rules.) 



3 The Unitarity Triangle 

The allowed region in p-r] space can be displayed quite elegantly using the so-called uni- 
tarity triangle. The unitarity of the CKM matrix leads to the following relation: 

(22) v ud v* b + v cd v; b + v td v t t = . 

Using the form of the CKM matrix in Eq. [I], this can be recast as 

(23) At+Aifr 1 - 

which is a triangle relation in the complex plane (i.e. p-t] space), illustrated in Fig. [l|. 
Thus, allowed values of p and r\ translate into allowed shapes of the unitarity triangle. 

In order to find the allowed unitarity triangles, the computer program MINUIT is used 
to fit the CKM parameters A, p and r\ to the experimental values of |^WK&|; l e l 
and Xd- Since A is very well measured, we have fixed it to its central value given above. 
As discussed in the introduction, we present here two types of fits: 

• Fit 1: the "experimental fit." Here, only the experimentally measured numbers are 



used as inputs to the fit with Gaussian errors; the coupling constants fB d yB Bd and 
Bk are given fixed values. 

Fit 2: the "combined fit." Here, both the experimental and theoretical numbers are 
used as inputs assuming Gaussian errors for the theoretical quantities. 
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We first discuss the "experimental fit" (Fit 1). The goal here is to restrict the allowed 
range of the parameters (p, rj) for given values of the coupling constants f Bd \J B Bd and 

Bk- For each value of Bk and f Bd \J ' B Bd , the CKM parameters A, p and r\ are fit to the 
experimental numbers given in Table p] and the x 2 is calculated. 

First, we fix B K = 0.8, and vary f Bd \f^B d in the range 130 MeV to 230 MeV. The 
fits are presented as an allowed region in p-rj space at 95% C.L. (x 2 = Xmin + 6-0). The 
results are shown in Fig. |2|. As we pass from Fig. 0(a) to Fig. @(e), the unitarity triangles 
represented by these graphs become more and more obtuse. Even more striking than this, 
however, is the fact that the range of possibilities for these triangles is quite large. There 
are two things to be learned from this. First, our knowledge of the unitarity triangle is at 
present rather poor. This will be seen even more clearly when we present the results of Fit 
2. Second, unless our knowledge of hadronic matrix elements improves considerably, mea- 
surements of |e| and Xd, no matter how precise, will not help much in further constraining 
the unitarity triangle. This is why measurements of CP- violating rate asymmetries in the 
B system are so important [f42| , fi3"| . Being largely independent of theoretical uncertainties, 
they will allow us to accurately pin down the unitarity triangle. With this knowledge, 

we could deduce the correct values of B K and f Bd \J 5^, and thus rule out or confirm 
different theoretical approaches to calculating these hadronic quantities. 

Despite the large allowed region in the p-r] plane, certain values of Bk and fB d \J B Bd 
are disfavoured since they do not provide a good fit to the data. For example, fixing 
Bk = 1-0, we can use the fitting program to provide the minimum x 2 f° r various values of 

f Bd \]~B~B~ d - The results are shown in Table ||, along with the best fit values of (p, if) . Since 
we have two variables (p and rj), we use Xmin < 2-0 as our "good fit" criterion, and we see 
that f Bd \j B Bd < 130 MeV and f Bd \f B Bd > 270 MeV give poor fits to the existing data. 
Note also that the x 2 distribution has two minima, at around f Bd \J 1 B Bd = 150 and 230 

MeV. We do not consider this terribly significant, since the surrounding values of f Bd yB Bd 
also yield good fits to the data. The very small values of x 2 m in depend sensitively on the 
central values of the various experimental quantities - if these values move around a little 
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Figure 2: Allowed region in p-rj space, from a fit to the experimental values given in Table 
[l]. We have fixed B K = 0.8 and vary the coupling constant product f Bd J ' B Bd as indicated 
on the figures. The solid line represents the region with x 2 — Xmin + 6 corresponding to 
the 95% C.L. region. The triangles show the best fit. 
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IP) 


V 2 

A mm 


120 


(-0.43, 0.13) 


2.89 


130 


(-0.38, 0.16) 


1.16 


140 


(-0.34, 0.18) 


0.3 


150 


(-0.30, 0.21) 


9.2 x 10~ 3 


160 


(-0.25, 0.24) 


0.07 


170 


(-0.20, 0.27) 


0.29 


180 


(-0.14, 0.29) 


0.52 


190 


(-0.07, 0.31) 


0.64 


200 


(0.0, 0.31) 


0.59 


210 


(0.07, 0.31) 


0.38 


220 


(0.13, 0.31) 


0.13 


230 


(0.19, 0.3) 


3.8 x 10~ 3 


9/i n 


fn oq n 3 1 ^ 

yU.ZOj U.Ol ) 


U.U ( 


250 


(0.27, 0.31) 


0.36 


260 


(0.31, 0.31) 


0.87 


270 


(0.34, 0.31) 


1.59 


280 


(0.38, 0.32) 


2.51 



Table 2: The "best values" of the CKM parameters (p,rj) as a function of the coupling 
constant fs d J Bs d , obtained by a minimum x 2 fit to the experimental data, including the 
renormalized value of m t = 170 ± 11 GeV. We fix B % = 1.0. The resulting minimum % 2 
values from the MINUIT fits are also given. 



bit, the values of fB d yE>B d which give the minimum x 2 values will move around as well. 
In Tables |] and |J we present similar analyses, but for Bk = 0.8 and 0.6, respectively. 
From these tables we see that the lower limit on fs d \J Bs d remains fairly constant, at 
around 130 MeV, but the upper limit depends quite strongly on the value of B K chosen. 



Specifically, for B K = 0.8 and 0.6, the maximum allowed value of fB d \jBB d is about 240 
and 210 MeV, respectively. 



In Table |5], we present the x 2 values as a function of fB d \j Bs d for Bk = 0.4, which is 
not favoured by lattice calculations or QCD sum rules. What is striking is that, over the 

entire range of fs d \f Bs d , the minimum x 2 is always greater than 2. This indicates that 
the data strongly disfavour Bk < 0.4 solutions. 

We now discuss the "combined fit" (Fit 2). Since the coupling constants are not 
known and the best we have are estimates given in the ranges in Eqs. (|16D and (|2l), 



a reasonable profile of the unitarity triangle at present can be obtained by letting the 
coupling constants vary in these ranges. The resulting CKM triangle region is shown in 
Fig. [| As is clear from this figure, the allowed region is rather large at present. The 
preferred values obtained from the "combined fit" are 

(24) (p, rj) = (-0.07, 0.34) (with X 2 = 6.6 x 10~ 2 ) . 
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fB^B Bd (MeV) 


(p,v) 


2 

X,min 


120 


(-0.42, 0.16) 


3.04 


130 


(-0.37, 0.19) 


1.32 


140 


(-0.32, 0.23) 


0.43 


150 


(-0.27, 0.26) 


0.07 


160 


(-0.22, 0.29) 


1.4 x 10" 3 


170 


(-0.15, 0.32) 


0.05 


180 


(-0.08, 0.34) 


0.09 


190 


(-0.01, 0.35) 


0.06 


200 


(0.06, 0.35) 


0.01 


210 


(0.13, 0.35) 


0.02 


220 


(0.18, 0.35) 


0.2 


230 


(0.23, 0.35) 


0.61 


240 


(0.28, 0.35) 


1.29 


250 


(0.32, 0.35) 


2.22 



Table 3: The "best values" of the CKM parameters (p, rf) as a function of the coupling 
constant fs d \J Bs d , obtained by a minimum \ 2 fit to the experimental data, including the 
renormalized value of m t = 170 ±11 GeV. We fix B K = 0.8. The resulting minimum \ 2 
values from the MINUIT fits are also given. 



f Bd {K d (MeV) 


(p,v) 


v 2 . 
Amm 


120 


(-0.4, 0.21) 


3.39 


130 


(-0.35, 0.25) 


1.7 


140 


(-0.29, 0.29) 


0.78 


150 


(-0.22, 0.33) 


0.35 


160 


(-0.15, 0.36) 


0.18 


170 


(-0.07, 0.38) 


0.16 


180 


(0.01, 0.39) 


0.24 


190 


(0.08, 0.4) 


0.48 


200 


(0.15, 0.4) 


0.96 


210 


(0.21, 0.4) 


1.73 


220 


(0.26, 0.4) 


2.85 



Table 4: The "best values" of the CKM parameters (p, rj) as a function of the coupling 

constant fs d \J Bb 4 , obtained by a minimum x 2 fit to the experimental data, including the 

renormalized value of mt = 170 ±11 GeV. We fix Bk = 0.6. The resulting minimum x 2 
values from the MINUIT fits are also given. 
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f Bd 4^ d (MeV) 


(p,v) 


v 2 . 

Amm 


130 


(-0.28, 0.35) 


2.95 


140 


(-0.2, 0.39) 


2.22 


150 


(-0.11, 0.43) 


2.04 


160 


(-0.02, 0.45) 


2.32 


170 


(0.06, 0.46) 


3.07 



Table 5: The "best values" of the CKM parameters (p, rj) as a function of the coupling 
constant fs d \[&~B d , obtained by a minimum \ 2 fit to the experimental data, including the 
renormalized value of mt — 170 ± 11 GeV. We fix = 0.4. The resulting minimum \ 2 
values from the MINUIT fits are also given. 




Figure 3: Allowed region in p-rj space, from a simultaneous fit to both the experimental 
and theoretical quantities given in Table [l|. The theoretical errors are treated as Gaussian 
for this fit. The solid line represents the region with \ 2 — Xmin + 6 corresponding to the 
95% C.L. region. The triangle shows the best fit. 
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Figure 4: Allowed region in p-rj space, from a simultaneous fit to both the experimental 
and theoretical quantities given in Table [I], except that we take B K = 0.6 ± 0.2. The 
theoretical errors are treated as Gaussian for this fit. The solid line represents the region 
with x 2 — Xmin + 6 corresponding to the 95% C.L. region. The triangle shows the best 
fit. 

For comparison, we also show the allowed region in the (p, rf) plane for the case in 
which Bk = 0.6 ± 0.2 [Eq. ([171) 1, which is more favoured by chiral perturbation theory 
and QCD sum rules. The CKM triangle region is shown in Fig. |j. Clearly, there is not 
much difference between this figure and Fig. |3|. The preferred values obtained from this 
fit are 

(25) (p, 77) = (-0.05, 0.37) (with X 2 = 0.1) . 



4 x s and the Unitarity Triangle 



Mixing in the B° s -B® system is quite similar to that in the B^-B® system. The B° s -B® box 
diagram is again dominated by t-quark exchange, and the mass difference between the 
mass eigenstates AM, is given by a formula analogous to that of Eq. (EO): 



(26) 



AM, = ^- 2 M 2 w M Bs (f B B Ba ) f) Bs y t f 2 (y t )\V*V t 



tb\ 



Using the fact that |V^,| = |Vt s | (Eq. [I]), it is clear that one of the sides of the unitarity 
triangle, |V^/AV^|, can be obtained from the ratio of AM d and AM S , 



(27) 



AM S _ Vb s M Bs [flBi 



A My 



VB d M Bd (f B B_ 



B, 



V t 



ts 



Vt 



td 



All dependence on the t-quark mass drops out, leaving the square of the ratio of CKM 
matrix elements, multiplied by a factor which reflects 5 , ^7(3)fl av0 ur breaking effects. The 
only real uncertainty in this factor is the ratio of hadronic matrix elements. Whether 
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or not be used to help constrain the unitarity triangle will depend crucially on 

the theoretical status of the ratio /§ B Bs j j\ B Bd . In what follows, we will take £ s = 



{fBs^B Ba )/(f Bd \jB Bd ) = (1.16 ± 0.1), consistent with both lattice-QCD || and QCD 
sum rules (40|. (The SU(3)-breaking factor in AM s /AM d is £ 2 .) 

The mass and lifetime of the B s meson have now been measured at LEP and Tevatron 
and their present values are M Bs = 5370.0 ± 2.0 MeV and r(B s ) = 1.58 ± 0.10 ps |2§. 
We expect the QCD correction factor fj Bs to be equal to its Bd counterpart, i.e. fj Ba = 
0.55. The main uncertainty in x s (or, equivalently, AM S ) is now f Bs B Bs . Using the 
determination of A given previously, r Ba = 1.58 ± 0.10 (ps) and Wit = 171 ± 11 GeV, we 
obtain 

AM S = (13.1 ±2.8) , f ^ BBB no (ps)- 1 , 
v ; (230 MeV) 2 Ky J 

f'B a BB s 
(230 MeV) 2 ' 

The choice f Bs \J B Bs = 230 MeV corresponds to the central value given by the lattice- 
QCD estimates, and with this our fits give x s ~ 20 as the preferred value in the SM. 



(28) x s = (20.7 ±4.5) 



Allowing the coefficient to vary by ±2a, and taking the central value for f Ba yB Bs , this 
gives 

11.7 < x s < 29.7 , 
(29) 7.5 (ps)- 1 < AM S < 18.7 (ps)' 1 . 

It is difficult to ascribe a confidence level to this range due to the dependence on the 
unknown coupling constant factor. All one can say is that the standard model predicts 
large values for x s , most of which are above the present experimental limit x s > 8.8 
(equivalently AM S > 6.1 (ps)- 1 ) |0|. 

An alternative estimate of AM S (or x s ) can also be obtained by using the relation in 
Eq. (|27|). Two quantities are required. First, we need the CKM ratio |V^/Vy. In Fig. [| 
we show the allowed values (at 95% C.L.) of the inverse of this ratio as a function of 



f Bd y B Bd , for B K = 0.8 ± 0.2. From this one gets 

V ts 



(30) 2.8 < 



V* 



Id 



< 7.6 



The second ingredient is the SU(3)-breaking factor which we take to be C, s — 1-16 ± 0.1, 
or 1.1 < £ 2 < 1.6. The result of the CKM fit can therefore be expressed as a 95% C.L. 
range: 

(31) 10.5 ( -V) < ^ < 77.7 ( is ^ 



1.16 y ~ AM d - V 1 - 16 / 

Again, it is difficult to assign a true confidence level to AM S / AMd due to the dependence 
on £ s . The large allowed range reflects our poor knowledge of the matrix element ratio 
| Vt s /Vtd| , which shows that this method is not particularly advantageous at present for 
the determination of the range for AM S . 
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Figure 5: Allowed values of the CKM matrix element ratio (V^/V^l as a function of the 

coupling constant product fs d \J ^B Bd , for B K = 0.8 ±0.2. The solid line corresponds to the 
best fit values and the dotted curves correspond to the maximum and minimum allowed 
values at 95 % C.L. 



The ALEPH lower bound AM S > 6.1 (ps)- 1 (95% C.L.) || and the present world av- 
erage AM,; = (0.465 ± 0.024) (ps)' 1 can be used to put a bound on the ratio AM S / AMd- 
The lower limit on AM S is correlated with the value of f s , the fraction of b quark fragment- 
ing into B s meson, as shown in the ALEPH analysis ||20|| . The value obtained from the 
measurement of the quantity f s BR(B s — > D s tvi) is f s = (12 ± 3)%. The time-integrated 
mixing ratios x an d Xd-, assuming maximal mixing in the B s -B s system Xs — 0.5, give 
f 8 = (9 ± 2)%. The weighted average of these numbers is f a = (10 ± 2)% @. With 
f s = 10%, one gets AM S > 5.6 (ps)- 1 at 95% C.L., yielding AMJ AM d > 11.3 at 95% 
C.L. Assuming, however, / s = 12% gives AM S > 6.1 (ps)' 1 , yielding AMJ AM d > 12.3 
at 95% C.L. We will use this latter number. 

The 95% confidence limit on AM S / AM d can be turned into a bound on the CKM 
parameter space (p,rf) by choosing a value for the SU(3)-breaking parameter £ 2 . We 
assume three representative values: £ 2 = 1.1, 1.35 and 1.6, and display the resulting 
constraints in Fig. ^. From this graph we see that the ALEPH bound marginally restricts 
the allowed p-rj region for small values of £ 2 , but does not provide any useful bounds for 
larger values. 

Summarizing the discussion on x s , we note that the lattice-QCD-inspired estimate 



fB s \lBB s — 230 MeV and the CKM fit predict that x s lies between 12 and 30, with a 
central value around 20. The upper and lower bounds and the central value scale as 



(/b s V-Bb s /230 MeV) 2 . The present constraints from the lower bound on x s on the CKM 
parameters are marginal but this would change with improved data. In particular, one 
expects to reach a sensitivity x s ~ 15 (or AM S ~ 10 ps -1 ) at LEP combining all data 



and tagging techniques [19|, which would be in the ball-park estimate for this quantity in 



the SM presented here. Of course, an actual measurement of x s would be very helpful in 
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Figure 6: Further constraints in p-rj space from the ALEPH bound on AM S . The bounds 
are presented for 3 choices of the SU(3)-breaking parameter: ^ = 1.1 (dotted line), 1.35 
(dashed line) and 1.6 (solid line). In all cases, the region to the left of the curve is ruled 
out. 



further constraining the CKM parameter space. 



5 CP Violation in the B System 

It is expected that the B system will exhibit large CP-violating effects, characterized 
by nonzero values of the angles a, (3 and 7 in the unitarity triangle (Fig. [1]) The 
most promising method to measure CP violation is to look for an asymmetry between 
T(B° — > f) and T(B° — > /), where / is a CP eigenstate. If only one weak amplitude 
contributes to the decay, the CKM phases can be extracted cleanly (i.e. with no hadronic 

uncertainties). Thus, sin 2a, sin 2(3 and sin 27 can in principle be measured in B d — > 7r + 7r~, 

Bd — > J/4>Ks an d B s — > pKg, respectively. 

Unfortunately, the situation is not that simple. In all of the above cases, in addition 



to the tree contribution, there is an additional amplitude due to penguin diagrams |44 . 
In general, this will introduce some hadronic uncertainty into an otherwise clean measure- 
ment of the CKM phases. In the case of Bd — > J/ipKs, the penguins do not cause any 
problems, since the weak phase of the penguin is the same as that of the tree contribution. 
Thus, the CP asymmetry in this decay still measures sin 2f3. 

For Bd — ► 7r + 7r~, however, although the penguin is expected to be small with respect 
to the tree diagram, it will still introduce a theoretical uncertainty into the extraction 



of a. Fortunately, this uncertainty can be removed by the use of isospin [|5]. The 
key observation is that the 1 = 2 component of the B — > titt amplitude is pure tree 
(i.e., it has no penguin contribution) and therefore has a well-defined CKM phase. By 
measuring the rates for B + — > 7r + 7r°, B° — > 7r + 7r~ and B° — > 7r°7r°, as well as their CP- 
conjugate counterparts, it is possible to isolate the 1 = 2 component and obtain a with 
no theoretical uncertainty. Thus, even in the presence of penguin diagrams, sin 2a can 
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Figure 7: Allowed region of the CKM matrix elements \Vtd\ and \V u b\ resulting from the 



'combined fit" of the data for the ranges for f Bd \lB Bd and B K given in the text. 



in principle be extracted from the decays B — > tttt. It must be admitted, however, that 
this isospin program is ambitious experimentally. If it cannot be carried out, the error 
induced on sin 2a is of order \P/T\, where P (T) represents the penguin (tree) diagram. 
The ratio \P/T\ is difficult to estimate - it is dominated by hadronic physics. However, 
one ingredient is the ratio of the CKM elements of the two contributions: |V^T4d/T^T4d| — 
l^td/Ktbl- From our fits, we have determined the allowed values of \V t d\ as a function of 
\V u b\. This is shown in Fig. |7| for the "combined fit". The allowed range for the ratio of 
these CKM matrix elements is 

14* I 



(32) 



1.2 < 



ub 



< 5.8 



with the central value close to 3. 

It is B s pKs which is most affected by penguins. In fact, the penguin contribution 
is probably larger in this process than the tree contribution. This decay is clearly not 
dominated by one weak (tree) amplitude, and thus cannot be used as a clean probe of 
the angle 7. Instead, two other methods have been devised, not involving CP-eigenstate 

final states. The CP asymmetry in the decay B s — > DfK T can be used to extract sin 2 7 
46fl . Similarly, the CP asymmetry in B^ — > D° CP K ± also measures sin 2 7 H47 |. Here, D° 



CP 



is a D° or D° which is identified in a CP-eigenstate mode (e.g. tt + it~, K + K~, ...). 

These CP-violating asymmetries can be expressed straightforwardly in terms of the 
CKM parameters p and r]. The 95% C.L. constraints on p and rj found previously can 
be used to predict the ranges of sin 2a, sin 2(3 and sin 2 7 allowed in the standard model. 
The allowed ranges which correspond to each of the figures in Fig. |2|, obtained from Fit 
1, are found in Table |[ In this table we have assumed that the angle j3 is measured in 

Bd — * J/^Ks, and have therefore included the extra minus sign due to the CP of the 
final state. 

Since the CP asymmetries all depend on p and rj, the ranges for sin 2a, sin 2(3 and 
sin 2 7 shown in Table ^| are correlated. That is, not all values in the ranges are allowed 
simultaneously. We illustrate this in Fig. Rl corresponding to the "experimental fit" (Fit 
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/ Bd ^(MeV) 


sin 2a 


sin 2(3 


sin 2 7 


1 

lOU 

155 
180 
205 
230 


n ap. n as 

— U.oo 

0.75 - 1.0 
-0.59 - 1.0 
-0.96 - 0.92 
-0.98 - 0.6 


u.zi — U.O 1 
0.31 - 0.56 
0.42 - 0.73 
0.49 - 0.86 
0.57-0.93 


u. iz — u.oy 
0.34 - 0.92 
0.68 - 1.0 
0.37- 1.0 
0.28 - 0.97 



Table 6: The allowed ranges for the CP asymmetries sin 2a, sin 2(3 and sin 2 7, corre- 
sponding to the constraints on p and rj shown in Fig. 0. Values of the coupling constant 

fs d \J Bs d are stated. We fix Bk = 0.8. The range for sin 2/3 includes an additional minus 
sign due to the CP of the final state J/^K S . 



1), by showing the region in sin2a-sin2/3 space allowed by the data, for various values of 



fB d \jB Bd . Given a value for f Bd \J B Bd , the CP asymmetries are fairly constrained. How- 
ever, since there is still considerable uncertainty in the values of the coupling constants, a 
more reliable profile of the CP asymmetries at present is given by our "combined fit" (Fit 
2), where we convolute the present theoretical and experimental values in their allowed 
ranges. The resulting correlation is shown in Fig. || From this figure one sees that the 
smallest value of sin 2/3 occurs in a small region of parameter space around sin 2a ~ 0.4- 
0.6. Excluding this small tail, one expects the CP-asymmetry in Bd — > J/^/Ks to be at 
least 30%. 

It may be difficult to extract 7 using the techniques described above. First, since 
B s —* DfK T involves the decay of B s mesons, such measurements must be done at 
hadron colliders. At present, it is still debatable whether this will be possible. Second, the 
method of using B ± — > D^, P K to obtain 7 requires measuring the rate for B + — > D°K + . 
This latter process has an expected branching ratio of < O(10 -6 ), so this too will be hard. 

Recently, a new method to measure 7 was proposed Using a flavour SU(3) 

symmetry, along with the neglect of exchange- and annihilation-type diagrams, it was 
suggested that 7 could be found by measuring rates for the decays B + — > n°K + , B + — > 
7T + K°, B + — > 7r + 7r°, and their charge-conjugate processes. The irK final states have both 
1=1/2 and 1 = 3/2 components. The crucial ingredient is that the gluon- mediated 
penguin diagram contributes only to the 1 = 1/2 final state. Thus, a linear combination 
of the B + — > 7T° K + and B + — > tt + K° amplitudes, corresponding to / = 3/2 in the irK 
system, can be related via flavour SU(3) to the purely 1 = 2 amplitude in B + — > 7r + 7r°, 
permitting the construction of an amplitude triangle. The difference in the phase of the 
B + — > 7T + 7r° side and that of the corresponding triangle for B~ decays was found to be 
27. SU(3) breaking can be taken into account by including a factor fx/fn in relating 
B — > 7T7T decays to the B — > ttK decays [f49f . 

The key assumption is that the penguin is mediated by gluon exchange. However, 
there are also electroweak contributions to the penguins |5(J. These electroweak penguins 
(EWP's) are not constrained to be isosinglets. Thus, in the presence of EWP's, there is 
no longer a triangle relation B — > ttK and B — > tttt amplitudes [[51]. Indeed, electroweak 
penguins can, in principle, even invalidate the isospin analysis in B — > tt-k, since the 1 = 2 
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Figure 8: Allowed region of the CP asymmetries sin 2a and sin 2/3 resulting from the 
"experimental fit" of the data for different values of the coupling constant f Bd \J B Bi 
indicated on the figures a) - e). We fix B K = 0.8. 
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= 180 + 50 MeV , B D = 1 .0 + 0.2 



= 0.8 + 0.2 



-0.4 -0.2 0.2 
sin 2a 



Figure 9: Allowed region of the CP asymmetries sin 2a and sin 2(3 resulting from the 



"combined fit" of the data for the ranges for fB d \jBs d and Bk given in the text, 
amplitude will include a contribution from EWP's, and hence will no longer have a well- 



defined weak CKM phase. However, theoretical estimates |51 , |52|] show that electroweak 
penguins are expected to be relatively unimportant for B — > tttt. 

The question of the size of EWP's has therefore become a rather interesting question, 
and a number of papers have recently appeared discussing this issue |53| . These include 



both theoretical predictions, as well as ways of isolating EWP's experimentally. The 
general consensus is that EWP's are large enough to invalidate the method of Ref. |48| for 



obtaining 7. However, two new methods making use of the flavour SU(3) symmetry, and 
which do not have any problems with electroweak penguins, have been suggested. Both 
are rather complicated, making use of the isospin quadrangle relation among B — > ttK 
decays, as well as B + — > tt + 7t° plus an additional decay: B s — > r]7c° in one case |52 |, 



B + — > r/K + in the other [5"3|. Although electroweak penguins do not cause problems, 
SU(3)-breaking effects which cannot be parametrized simply as a ratio of decay constants 
are likely to introduce errors of about 25% into both methods. It is clear that this is a 
subject of great interest at the moment, and work will no doubt continue. 



6 Summary and Outlook 

We summarize our results: 

(i) We have presented an update of the CKM unitarity triangle using the theoretical 
and experimental improvements in the following quantities: \V c b\, \V u b/V c b\, AM d , r(B d ), 
Tnt, Vcc Vet- The fits can be used to exclude extreme values of the pseudoscalar coupling 
constants, with the range 130 MeV < f Bd \[B~^ d < 270 MeV still allowed for B K = 1. The 
lower limit of this range is quite -Bjf-independent, but the upper limit is strongly correlated 
with the value chosen for B K . For example, for B K = 0.8 and 0.6, fs d \J Bs d < 240 and 

210 MeV, respectively, is required for a good fit. The solutions for Bk = 0.8 ± 0.2 
are slightly favoured by the data as compared to the lower values. These numbers are 
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in very comfortable agreement with QCD-based estimates from sum rules and lattice 
techniques. The statistical significance of the fit is, however, not good enough to determine 
the coupling constant more precisely. We note that B^ < 0.4 is strongly disfavoured by 
the data, since the quality of fit for such values is very poor. 

(ii) The newest experimental and theoretical numbers restrict the allowed CKM uni- 
tarity triangle in the (p, ?7)-space somewhat more than before. However, the present 
uncertainties are still enormous - despite the new, more accurate experimental data, our 
knowledge of the unitarity triangle is still poor. This underscores the importance of mea- 
suring CP-violating rate asymmetries in the B system. Such asymmetries are largely 
independent of theoretical hadronic uncertainties, so that their measurement will allow 
us to accurately pin down the parameters of the CKM matrix. Furthermore, unless our 
knowledge of the pseudoscalar coupling constants improves considerably, better measure- 
ments of such quantities as Xd will not help much in constraining the unitarity triangle. 
On this point, help may come from the experimental front. It may be possible to measure 
the parameter fs d , using isospin symmetry, via the charged-current decay —>■ r ± v T . 
With \V ub /V cb \ = 0.08 ± 0.02 and f Bd = 180 ± 50 MeV, one gets a branching ratio 
BR(B± r ± z/ T ) = (1.5-14.0) x 10~ 5 , with a central value of 5.2 x 10~ 5 . This lies in the 
range of the future LEP and asymmetric I?-factory experiments, though at LEP the rate 
Z — > B C X — > t u T X could be just as large as Z — > B^X — > t u T X. Along the same 
lines, the prospects for measuring (Jb^/bJ m the FCNC leptonic and photonic decays 
of B® and B® hadrons, (B%, B®) — > //+ pT , [B% B®) — > 77 in future B physics facilities are 
not entirely dismal |55|j . 

(iii) We have determined bounds on the ratio |Vtd/Vf S | from our fits. For 130 MeV < 



f Bd y B Bd < 270 MeV, i.e. in the entire allowed domain, at 95 % C.L. we find 

V td 



(33) 0.13 < 



Vt 



tH 



< 0.35 



The upper bound from our analysis is more restrictive than the current experimental upper 
limit following from the CKM-suppressed radiative penguin decays BR(B — > uj + 7) and 



BR(B — > p + 7), which at present yield at 90% C.L. |56 

V td 



(34) 



V t 



ts 



< 0.64-0.75 



depending on the model used for the SU(3)-breaking in the relevant form factors [57 . 
Long-distance effects in the decay B^ — ► p ± + 7 may introduce theoretical uncertainties 
comparable to those in the SU(3)-breaking part but the corresponding effects in the decays 
B° — > (p°,uj) + 7 are expected to be very small [f)8|| . Furthermore, the upper bound is 



now as good as that obtained from unitarity, which gives 0.08 < (V^/Vy < 0.36, but the 
lower bound from our fit is more restrictive. 

(iv) Using the measured value of mt, we find 

f 2 B 

(35) x s = (20.7 ± 4.5) ' /! H! ~ 



(230 MeV) 2 
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Taking fB a \jBB s = 230 (the central value of lattice-QCD estimates), and allowing the 
coefficient to vary by ±2o~, this gives 

(36) 11.7 < x s < 29.7 . 

No reliable confidence level can be assigned to this range - all that one can conclude is 
that the SM predicts large values for x s , most of which lie above the ALEPH 95% C.L. 
lower limit of x s > 8.8. 

(v) The ranges for the CP-violating rate asymmetries parametrized by sin 2a, sin 2/5 
and and sin 2 7 are determined at 95% C.L. to be 

-1.0 < sin 2a < 1.0 , 

(37) 0.21 < sin 2(3 < 0.93 , 

0.12 < sin 2 7 < 1.0 . 

(For sin 2a < 0.4, we find sin 2/3 > 0.3.) Electroweak penguins may play a significant role 
in some methods of extracting 7. Their magnitude, relative to the tree contribution, is 
therefore of some importance. One factor in determining this relative size is the ratio of 
CKM matrix elements | V*rf / V^fe | . We find 



(38) 1.2 < -f- < 5. 



ub 
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